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Abstract 

A fully anisotropic simple-cubic Ising lattice in the geometry of periodic cylin- 
ders nxnxoois investigated by the transfer-matrix finite-size scaling method. 
In addition to the previously obtained critical amplitudes of the inverse cor- 
relation lengths and singular part of the free energy [M. A. Yurishchev, Phys. 
Rev. B 50, 13 533 (1994)], the amplitudes of the usual ("linear") and nonlin- 
ear susceptibilities and the amplitude of the second derivative of the spin-spin 
inverse correlation length with respect to the external field are calculated. The 
behavior of critical amplitude combinations (which, in accordance with the 
Privman-Fisher equations, do not contain in their composition the nonuniver- 
sal metric coefficients and geometry prefactor) are studied as a function of the 
interaction anisotropy parameters. A universality domain for the amplitude 
ratios is found in the quasi-one-dimensional regime of interactions in the sys- 
tem. In the case of a fully isotropic three-dimensional Ising model for which 
the high precision values of the critical coupling and critical-point free energy 
are available, improved estimates are obtained for the following four universal 
quantities: (1) the amplitude of spin-spin inverse correlation length, (2) the 
amplitude of singular part of the free energy, (3) the ratio of the amplitude 
of a second derivative of the spin-spin inverse correlation length with respect 



to the external field to the usual susceptibility amplitude, and (4) the ratio of 
the nonlinear susceptibility amplitude to the square of the linear susceptibility 
amplitude (i.e., for the finite-size counterpart of the four-point renormalized 
coupling constant). 

PACS numbers: 05.50.+q, 05.70.Jk, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 



In paper JIJ (hereafter referred to as I), the critical finite-size amplitudes of the inverse 
correlation lengths and the "singular" part of the free energy per spin have been calculated 
for the three-dimensional ferromagnetic Ising model on a simple-cubic lattice with fully 
anisotropic interactions. The amplitudes have been obtained by a finite-size scaling (FSS) 
analysis of the transfer matrix (TM) data for subsystems having a shape of n x n x oo 
cylinders with periodic boundary conditions in both transverse directions. 

The results showed that the ratios of the above critical amplitudes practically do not 
depend on the anisotropy parameter J x j J z if the second parameter J y / J x is fixed (the 
parallelepipeds n x n x oo are assumed to be stretched along z axis; J x and J y are the 
interaction constants in the transverse directions of a cylinder, and J z is the interaction 
constant in the longitudinal direction). Furthermore, it has been found that the amplitude 
ratios are independent of the second anisotropy parameter, J y / J x , in the vicinity of J y / J x = 
1. Unfortunately, a wrong conclusion concerning the amplitude ratio constancy in the full 
range of the J y / J x values has been made in I. 

In addition to the amplitude of the spin-spin inverse correlation length (A s ) and the 
free-energy amplitude (Af), in the present paper we calculate the amplitudes of linear and 
nonlinear susceptibilities (A x and A x (±), respectively), as well as the amplitude of the second 
derivative of the spin-spin inverse correlation length with respect to the external field (A K ^ h ). 
The calculations are performed again for the cyclic clusters n x n x oo with the maximum 
number of lattice layers in each transverse direction n = 4 (i.e., for the sizes of TMs up to 
65 536 x 65 536). Guiding ourselves by the functional equations of Privman and Fisher . 
we construct from critical amplitudes the combinations A K » / ' A x and A x {a)A s /A x which, as 
A s /Af, do not contain the certain nonuniversal coefficients. The analysis shows that those 
amplitude ratios are universal with respect to the anisotropy parameter J x / J z — > but 
depend on the parameter J y /J x . However, the J y /J x dependencies exhibit smooth extrema 
near J y j J x — 1, and therefore the amplitude ratios have a local universality in their vicinity 
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upon the second anisotropy parameter of the model. 

In the present paper we find also the values of A s , Af, A K " /A x , and j A 2 X for n — 2, 3, 
and 4 in the special case of the fully isotropic simple-cubic Ising lattice by using the available 
high accurate estimates for the critical point and critical free energy and give three-point 
extrapolations for these quantities. 

II. SPATIAL ANISOTROPY AND THE PRIVMAN-FISHER EQUATIONS 

In the case of d- dimensional cylinders with sizes of L\ x L 2 x ■ ■ ■ x L^i x oo, the 
Privman-Fisher equations which allow to identify for the hyperscaling systems the universal 
combinations of critical finite-size amplitudes have the form || 

k(t, h) = L^X{C x tLl T , C 2 hLf- LjLo, L d ^/L ) (1) 

and 

/«(*, h) = L d Y(C±tLl T , C 2 hLl h ;Li/Lo, L d ^/L ) . (2) 

In these equations, k is the inverse correlation length in the longitudinal direction of the 
cylinder, f( s > is the singular part of the free-energy density measured in units of — /cgT, 
t — (T — T c )/T c is the reduced temperature, h is the normalized external field, dt and 
are critical exponents, C\ and C 2 are system-dependent metric factors, and Lq is a scale 
length. The FSS functions X(x-\_, x 2 \ h, l 2 , . . .) and Y(x\, x 2 \ l\, l 2 , . . .) are universal (within 
the limits of a given universality class) but may depend on the type of boundary conditions 
and, as we see, on the aspect ratios k = L^/Lq. For the scaling length L , one can take 
L Q = {L\Li . . . or put L = min(Li, L 2 , . . . , etc. 

Let us pass in Eqs. and (|2|) to the dimensionless quantities which are immediately 
defined by the parameters of a subsystem Hamiltonian. Let the elementary cell of a simple- 
cubic lattice have the sizes of a\ x a 2 x ■ ■ • x ad and let the subsystem "frame" be a cylinder 
rix x n 2 x ■ ■ ■ x rid-i x oo, where is a number of spins in the zth transverse direction. 
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We will, however, restrict ourselves to the case of cylinders with square cross section, i.e., 
m — n 2 — ... = n d -i = n. Then Lj = na,-, « = a^" 1 /^, and = (oi0 2 • ■ ■ a d )~ l f^\ 
where n n and are, respectively, the dimensionless inverse correlation length and singular 
part of the free energy per site. (Below, however, we will often omit, for brevity, the word 
"dimensionless".) In the new variables the Privman-Fisher equations can be written as 

K n (t, h) = n - 1 GX(C 1 tn VT , C 2 hn Vh ;a 2 /ai, . . . , a rf _i/ai) (3) 

and 

/( s )(i,/ 1 ) = rtr(C 1 in»C^;a 2 /a ll ...,a H /a 1 ) , (4) 

where G = G(a 2 /ai, . . . , a d _i/ai, a d /ai) is the geometry prefactor. We suppose that at least 
two lattice spacings in the system are finite and nonzero; without loss of generality, we take 
them as a\ and a d . 

One obtains from Eqs. (|3|) and (^) that 

A s = nn n (0, 0) = GX(0, 0; a 2 /a u . . . , a d -i/ cti) , (5) 

Af = n d fi s \0, 0) = G F(0, 0; a 2 /a h . . . , a^/a^ , (6) 
and analogously for the dimensionless critical finite-size-scaling amplitudes of the derivatives 

Of K n (t,h) Ovf( s \t,h). 

The geometry-pref actor form depends on the choice of L . The prefactor is defined 
with exactness up to some multiplicative function of a 2 /a\, . . . ,a d _i/a\, inasmuch as such 
a function can be introduced into the scaling functions X and Y or, vice versa, taken out 
from them. We will however put G = 1 for the fully isotropic model. 

In order to be able to employ the FSS theory, the "lattice spacings" Oj must be chosen so 
that in the vicinity of the phase transition point the bulk correlation lengths £i, • • • , £d along 
all the different spatial directions become equal between themselves: £i = £2 = • • • = £d- 
Such a rescaling of lattice spacings can be made for systems with the isotropic critical expo- 
nent v of the bulk correlation lengths, i.e. when = ($t~ v (Refs. |5|-0). The choice of ai 
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completes, in principle, the process of expressing the quantities entering into the Privman- 
Fisher equations via "microscopic" parameters — the interaction constants of the Hamil- 
tonian. Note that the situation is paradoxical to a certain extent: in order to use the FSS 
theory for extracting information about a bulk system from properties of its finite subsys- 
tems, one must know before the correlation-length amplitudes of the bulk system. 

In the two-dimensional anisotropic Ising lattice which is a limited case (J y = 0) of the 
three-dimensional Ising model under question, the dimensionless bulk correlation lengths for 
T >T C are § (see also, e.g., Ref. [|) 

U-Jx/k B T, J z /k B T) = U-Jz/k B T, J x /k B T) = 1/ hx[coth{J x /k B T) exp(-2 J z /k B T)] . (7) 

(When T < T c , the correlation length expressions differ from Eq. (|7|) by the factor —1/2 
which is, however, unessential for the present considerations.) According to the isotropy 
requirement, a x £ x = a z ^ z and therefore for the square Ising lattice [f| (see also [|6|, |lO| , pT| ) 



G 



sinh 



2J Z 
k B T c 



(8) 



where the critical temperature T c is determined by the equation 

2 J f \ . / 2 J ■ 
smh 



sinh 



1 . 



(9) 



k B T c J \k B T c 

Eqs. (H|) and (^J) in parametric form yield the geometry factor G(J X /J Z ). Note that up to 
date the anisotropy factors have been found for many exactly solved two-dimensional models 
of statistical physics [P^H 



In the three-dimensional space, Eqs. (^) and (Q) are written as: 

K n (t, h) = n- 1 GX(C 1 tn VT ,C 2 hn yh ]a y /a x ) 



(10) 



and 



fi s) (t, h) = n-'- i GY{C l tn yT ,C 2 hn yh ]a y /a x ) . 



(11) 



Due to the shape parameter a y /a x , the scaling functions depend on the coupling ratios, 



and therefore Eqs. (|10|) and (|TT|) do not allow the existence of any universal amplitude 



combinations. An obvious exception, however, is given by the case when the interactions in 
the transverse directions of the parallelepiped n x n x oo are equal between themselves. On 
physical grounds, it is clear that a x = a y by this, and consequently the critical amplitude 
combinations which do not contain C\, C2, and G will not depend upon the interaction 
constants in the system. 

III. FORMULAS FOR THE SUSCEPTIBILITIES AND FOR THE INVERSE 
CORRELATION-LENGTH DERIVATIVE 

For the amplitude ratio A » /A x , it follows from Eqs. (|3]) and (|J) that 

hh 

A >, 



K hh 



(12) 



A x n d - l Xn 

where K hh = d 2 K n /dh 2 , while Xn = d 2 f n /dh 2 = d 2 f^/dh 2 is the susceptibility of the 
system (f n denotes the dimensionless free energy per lattice site). Taking into account that 
A s = nn n , one finds for the second amplitude combination 

A xW A s x { n ] Kn , v 

A\ n^xl ' 1 6) 

where Xn" 1 = —d^fn/dh 4 is a nonlinear susceptibility. Thus, the problem is to carry out the 
calculation of K hh and x^n f° r periodic cylinders nxnxooat/i = since the method of 
calculating n n and Xn w as already done in I. 

To solve these problems, we use, as in I, the TM technique. The matrix elements of TM, 
V, are given by 

n 

(Sn, S12, . . . , S nn \V\S' u , S' 12 , . . . , S' nn ) = Y\_ ex p[2-Kx(SijS i+ ij + S[jS' i+l j) 

+\K y (S ij S ij+ i + S'^S'^) + KzSijS'ij + \h(Sij + S'ij)] . (14) 

Here the spin variables Sij take the values ±1; Si n+ i = Sn and S n+ ij = Sy for all i, j = 
1,2, ... ,n; K a = Ja/ksT (a = x, y, z). The matrix V is real, symmetric, and dense; having 
all its elements positive. The dimensionless free energy per spin equals 



rr 



i • 



(15) 



where Ai is the largest eigenvalue of V. Further, the dimensionless spin-spin inverse corre- 
lation length in the longitudinal direction of a parallelepiped n x n x oo is 



n n = bi(A l /A 2 ) 



(16) 



where A2 denotes the second largest eigenvalue of the matrix V. 

In order to derive the exact formulas for the zero-field derivatives of the free energy and 
inverse correlation length, we will use perturbation theory. For this one expands the TM in 
powers of h: 



V = V + hV 1 + h 2 V 2 + h 3 V 3 + h*V 4 + 0{h 5 ) . 



(17) 



Let us use the symmetry of the model (in this context see I). V is invariant under the 
transformations of the group TAC 2v (T is a group of translations in the transverse directions 
of a cyclic bar n x n x 00 and C 2v is the point group generated by two symmetry planes 
going through the middles of opposite faces of the system). Pass in Eq. ( |TTD to the basis 
of the identity irreducible representation of the group T A C 2v . Expansion ([TtD preserves 
the above form but now its terms are blocks corresponding to the indicated representation. 
Both eigenvalues Ai and A2 lie in the given block V. Take now into consideration the 
symmetry Z 2 (a group of spin inversions). The matrices V, V 2 , and V4 are symmetrical and 
the matrices V\ and V3 are antisymmetrical (i.e., they change sign) under the spin inversion 
operation. Going by means of a similarity transformation into a new basis in which the 
original representation of the group is completely reducible, one obtains from Eq. (IT] 
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+h 4 
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v 4 



0(h 5 



(18) 
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For definiteness we suppose that the subblocks V 2 \ and V4 correspond to the identity 
irreducible representation of the group Z 2 . Denote the sizes of these subblocks by N± x N\ 
and the sizes of the subblocks V^ 2 \ V 2 , and v} 2 ^ by N 2 x A^ 2 - Then, the matrices V"/ 12 '* and 
V-} 12 ^ have the sizes of JVi x ]\f 2 ; v} 21 ^ and V^ 21 "* are found by transposing V"/ 12 "* and V^ 12 \ 
respectively. The numbers Ni and N 2 can be obtained from a group-theoretical analysis. 
For the 2 x 2 x 00 cluster, N x = 5 and N 2 = 2 (Ref. @ ) ; for the 3 x 3 x 00 one, N t = N 2 = 18 
(Ref. [T§); and for 4 x 4 x 00, JVi = 787 and N 2 = 672 (Ref. I). Let, further, A, a) and ^ be 
the eigenvalues and corresponding eigenvectors of matrix V^; let Aj = Ai be the largest 
eigenvalue of the block and consequently of the "nonperturbed" transfer matrix V. Due 
to the Perron theorem, Ai is nondegenerate. Finally, let \\ and tfi be the eigenpair of V^; 
Ai is the largest eigenvalue of block V^ 2 \ Note that A^ is also nondegenerate if we do not 
take the extreme cases which can be examined separately. 

Using the stationary perturbation theory for a nondegenerate level, we find the largest 
eigenvalue of V with accuracy up to the terms of second order in h: 



Ai = Ai + 



^ + vfVi + E 



k=l 



Ai - A*. 2 



/i 2 + 0(/i 4 ) . 



(19) 



From here, the expression for the initial (zero-field) susceptibility follows as 



Xn 



n 2 Xi 



^1^2 Vi + E 



fc=i 



Ai - A fe 2 



(20) 



The given formula, obtained by the perturbation-theory method, differs in a form from the 
one derived in I using the fluctuation-dissipation relation, but it is equivalent to it and yields 
the same values of the susceptibility. 

In an analogous way, one obtains the following expression for the second largest eigenvalue 
of V: 



A 2 



Ai 2) + 



\(2) 
k=l Ai 



A 



(1) 



h 2 + 0{h A ) . 



(21) 



From Eqs. flT5|), (JT9|), and (^Tj) we get the work formula for the second derivative of the 
spin-spin inverse correlation length at point h = 0: 



K 



hh 



Ai 



it=i Ai — A 



(2) 



Ai' 



(2) 



tftfV + E 



£ (^i (1 V) s 



^=1 



(2) 



(22) 



Calculating the largest eigenvalue Ai up to the terms of fourth order in h, we find the 
following result (suited for programming) for the initial nonlinear susceptibility: 



X 



(-1) 



12 



^Q 2 - 2(Q 1 + Q 2 + Q 3 + Q i + Q 5 -Q 6 + Q 7 - Q 8 ) 



(23a) 



with 



q = ^ + v 2 (1 Vi + E 



k=l M — A fc 



Qi = V^4 (1 Vi , Q2 = E 



^2 f^+JfWj, \2 



WTO 



fc=2 Ai — A fc 



ft = 2 fi^)(^) 



fc=l 



Ai - a[. 2) 



Q 4 = 2 E E ^ (12 V fc )(^ + ^ (21 V0(^ + ^ (1 Vi) 



fc=l Z=2 



(Ai - Ai 2) )(Ai - Aj iJ ) 



n _ ^ (^K (i2 y fc )(^ 2 (2 Vo(^ (2i Vi) 

"'"^ (Ai - Af } )(Ai - A[ 2) ) 



k=l 1=1 



q 6 = Vi v£ Vi E 



fe =i (A!-Ai 2) )2 



(23b) 



n _ f f f (^^n (12 V fc )(^^ 21 Vo(^^ 12 V m )(^ (21 Vi: 

7 (Ax-A^CAx-AP^Ax-Al?) 



k=l 1=2 m=l 



k=l [M — A fe ) i=i Ai — A; 

Explicit expressions for the matrix elements of subblocks vl l2 \ V-j; , V^, V^ 12 \ and 
are given in the Appendix. Notice that the full eigenproblems for the matrices and 



were solved by means of the C library function pair tred2 and tqli |16| . All calculations 
were carried out on a personal computer IBM PC-486 with the operating system FreeBSD. 
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IV. BEHAVIOR OF THE CRITICAL AMPLITUDE RATIOS 



As already noted in Sec. |TT], the FSS amplitude combinations normally should depend 
on the interaction ratios which change the shape of a subsystem. This is in agreement with 
basic ideas of the renormalization-group theory. 

We will consider how the amplitude combinations behave versus the interaction 
anisotropy parameters in the three-dimensional Ising model. To eliminate as effectively 
as possible nonuniversal quantities from the critical FSS amplitude combinations, it is natu- 
rally to consider the behavior of such combinations which do not include in its composition 
the nonuniversal metric and geometry factors. In other words, it is reasonable to base again 
on the Privman-Fisher equations in choosing the combinations. 

For a cluster nxnxoo, the values of linear and nonlinear susceptibilities and derivative 
K hh are taken at critical temperatures which were determined by the (n — l,n) cluster pair, 
i.e., through the equation 

(n - lK^Tj™- 1 ^) = n Kn (7f . (24) 

Since we are able to perform the calculations for subsystems with n < 4, this allows us 
to have two independent iterations, with (2,3) and (3,4) pairs (a degenerate pair (1,2) 
was eliminated from the consideration due to its anomalies [0]). Two steps of converging 
iterations give already a possibility to extract valuable information about tendencies in the 
change of quantities with increasing n. 

The results for the critical temperatures and critical FSS amplitude ratios A K " j A x and 
A x (^A s j ' A 2 X depending on J x / J z and J y / J x in the cases of both cluster pairs (2, 3) and (3, 4) 
are collected in Table I. 



Let us first discuss the behavior of amplitude ratios against J x /J z . Fig. 1 (T^j shows 
the dependencies of A » / A x on J x /J z at different fixed values of J y /J x . If J y /J x = 0, the 

hh 

bar n x n x oo is reduced to statistically independent strips. In two-dimensional space, the 
dependence on an anisotropy parameter must be, as it follows from Eqs. (Bf) and (Hf), absent. 
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The J y / J x = plot given in Fig. 1 confirms this: in the case of (3, 4) approximation, the 
ratio A n I A x is unchanged within the relative root-mean-square error 0.08%. We expect 

hh 

the constancy of amplitude ratios also at J y j J x — 1. The corresponding line shown in Fig. 1 
can be considered as a straight one with an accuracy 0.31%. Importantly, the analogous 
error for the (2,3) pair is 0.45%. Thus, the deviations versus J x /J z fall as n increases. 

For intermediate fixed values of J y /J x , the amplitude ratio has weak but still quite 
detectible dependence against J x /J z . As it is seen in Fig. 1, variances are most essential 
in the range 10 _1 ^ J y /J x < 1. Now the amplitude ratio dependence does not tend to 
disappear with increasing n. For example, when J y / J x = 0.1, the mean values and the root- 
mean-square uncertainties (the latter are shown in parentheses) for the quantity A » j A x 

hh 

in the range 10~ 3 < J x /J z < 1 equal 1.37(2) and 1.35(2) for the pairs (2,3) and (3,4), 
respectively. 

Similar picture is observed for the A X (^A S /A X (see Table I) and for the A s /Af (Table I 
of Ref. I). As a whole one can conclude that the critical FSS amplitude combinations which 
do not contain the nonuniversal factors C%, C2, and G display a tendency to the universality 
under the anisotropy parameter J x /J z , when this is small (i.e., by a quasi-one-dimensional 
nature of interactions in the system). 

We discuss now the behavior of the amplitude ratios as a function of the second anisotropy 
parameter, namely J y / J x . In I the amplitudes of the inverse correlation lengths and the free 
energy have been calculated for J y / J x G [0, 1]. The amplitude Af was found from a set of 
two equations 

f n = J 00 + n~ d 'A f (25) 

with n = 3 and 4 ["background" foo is a second unknown variable in Eq. (p5])]. In this paper 
we prolong such calculations to J y / J x = 4 and use not only the (3, 4) pair but also the (2, 3) 
one, to observe the evolution with increasing n. The results are shown in Fig. 2. It can be 
seen from the figure that when the parameter J y /J x increases from zero, the ratio A s /Af 
first increases monotonically from some finite value which decreases with growth of cluster 
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size. Then the A s /Af attains a smooth maximum near J y /J x = 1. Lastly, the ratio A s /Af 
falls monotonically as J y / J x becomes large. 

It follows from the obtained data that the values of A s /Af are equal between themselves 
for Jy/J x and (Jy/Jx)' 1 and this is better the smaller J x /J z . Clearly, the source of the 
Jy/Jx <-> (Jy/Jx)' 1 symmetry is caused by the (approximate) J x /J z independency of ampli- 
tude ratios. The inversion symmetry leads in turn to an existence of extremum at J y /J x = 1- 
Inasmuch as the extremum (maximum) is smooth, there is a local universality under the 
second anisotropy parameter J y /J x . We conclude that in the domain < J x /J z "C 1 and 
\Jyl Jx — 1| "C 1 there exist a complete universality of the amplitude ratio A s /Af. Note also 
that in the maximum region both curves go quite near one to the other; such neighborhood 
characterizes the arrived convergence of (2,3) and (3,4) approximations. 

We make the quantitative comparison of our data in the two-dimensional case. Space 
dimensionality d changes discontinuously from 3 to 2 in the limit J y /J x — > 0. This leads in 
turn to a finite jump in values of A s /Af at J y /J x = + and J y /J x = 0. Using Eq. (|25| ) one 
finds for the (n — 1, n) pair that 

/ „ i. s n(n — 1) (2n — 1) , . , , s 

(-4.M/>» = L(n - 1) + 1 ■ (*) 

where the subscripts and + denote the values at J y /J x = and + , respectively. Hence, 
recalculating the values at J y /J x = + to the values at J y /J x = 0, we obtain (A s /Af) = 
3.0143 for the (2,3) pair in the case J x /J z = 10~ 3 . For the (3,4) pair at the same value of 
J x /J z , (A s /Af) = 3.0087. In the two-dimensional Ising lattice, the exact value A s /Af = 3 
(in Fig. 2, it is shown by the open circle on the ordinate axis). Consequently, the error is 
reduced from 0.48% to 0.29% by going from (2, 3) to (3, 4) pair. 

Qualitatively similar picture takes place for the ratios A K " / A x and A X (4,)A S /A 2 X (Figs. 3 
and 4). Here again the amplitude ratios have finite values at J y /J x = + which vanish with 
increasing n according to the 1/n law. There is again a monotonic increase for small values 
of Jy/ J x , a broad maximum near J y / J x = 1 and then a monotonic decrease by large J y /J x . 
The inversion symmetry J y /J x <-» (Jy/Jx)' 1 applies again. 
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Let us compare the results for the discussed amplitude ratios in the two-dimensional 
limit. Since 

iK" / A x)o = n{A„ /A x ) 0+ , (27) 

lin nn 

taking from Table I the values at J y /J x = + , one finds that (A K >t /A x ) equals 1.959 and 

hh 

1.962 for the (2,3) and (3,4) approximations, respectively; (the point J x /J z = 1CT 3 has 
been used). These values are in good agreement with available estimates of the amplitude 
ratio for the two-dimensional isotropic Ising models, A K » /A x = 1.95 — 1.96 (Ref. ||19|| ). 

hh 

One can argue analogously for the combination A X (4)A S /A X , calculation being carried 
out by the formula 

(A xW A s /A 2 x ) = n(A xW A s /A 2 x ) 0+ . (28) 

Taking from Table I the necessary data, we find that at J x /J z = 10 -3 the ratio 
(A X (4)A S / A x ) = 5.7892 and 5.7976 for the pairs (2,3) and (3,4), respectively. These re- 
sults are in excellent agreement with the value A {a)A s /A 2 = 5.79728(5) following from the 



calculations IF20 



We discuss now a reason which could lead to a disappearance of the parameter J x /J z 
from scaling functions for the system with the dominant intrachain interaction (J z ^> J x , J y ). 
The bulk correlation-length amplitude is a function of the interaction constants. Let 
£o = <p(J x /kBT c , Jy/ksTc, J z /kBT c ). In a lattice which is infinite in all three directions, 
the function 0(x, y, z) is symmetrical under the replacement of its 2nd and 3rd arguments: 
4>(x,y,z) = <f)(x,z,y). Moreover, (p(0,y,z) = which means the absence of correlations in 
the direction along which there are no interactions. Let the expansion on the first argument 
begin from a term x in some power q. In accordance with physical reasoning, it is apparent 
that in the infinite lattice the amplitude ^ has the same functional form as but, of 
course, with replaced arguments: £q = 4>(J y /k B T C} J x /k B T C} J z /ksT c ). Further, the critical 
temperature of a quasi-one-dimensional Ising model is given by pif 



k B T c / J z = 2 



In ( -y^-r) ~ In In ( y^y) + 0(1) 



(29) 
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Due to the fact that a decrease of k B T c / J z is logarithmically slow, the arguments of function 
J x /k B T c = J x / J z ■ J z /k B T c and J y /k B T c = J y / J x ■ J x / J z ■ J z /k B T c are small for small J x / J z . 
Expanding the bulk correlation-length amplitudes in Taylor series and keeping only the 
leading asymptotical term, we obtain that as J X J J z — > the aspect ratio a y /a x ~ (J y / J x )~ q . 
The anisotropy parameter J x /J z did drop out. 

In the light of the above statements the Privman- Fisher equations for the quasi-one- 
dimensional system ( J x , J y -C J x ) can be written as 

Kn(t, h) = n- l G{J x /J z , J y /J x )X(C 1 tn yT ,C 2 hn y »; J y /J x ) (30) 

and 

fi s \t,h) = n- d G(J x /J z ,J y /J x )Y(C 1 tn yT ,C 2 hn^;JjJ x ) , (31) 

where d = |sign(J x )| + |sign(J 2/ )| + |sign(J 2 )|. The geometry pref actor is normalized so that 
G(l, 1) = 1 and the dependence G(J X / J z , 0) is given by Eqs. (§) and (|]). 

We know that the J y /J x <-> (J y / Jx)' 1 invariance and the analyticity of scaling functions 
leads to the existence of a smooth extremum for the critical amplitude ratios of the nxnxoo 
parallelepipeds at J y /J x = 1. In this context note that in the case of parallelepipeds with a 
rectangular cross section, n x X n y X oo, one should expect the extrema for the certain ratios 
of critical amplitudes at J y /J x = [n y jn, x ) x l <l . 



V. COMPLETELY ISOTROPIC LATTICE 

At present, the critical point of the fully isotropic simple-cubic Ising lattice is located 
to a high degree of accuracy: K c = 0.221655(1) (Ref. p2| and references therein). It is 



known also with large accuracy the free energy at criticality: = 0.77790(2), Ref []23 
On the other hand, we can carry out the calculations for subsystems n x n x oo with three 
numbers of lattice layers n — 2, 3, and 4. This allows, generally speaking, to perform the 
three-point extrapolations for accelerating the convergence and to improve the estimates of 
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the universal critical amplitudes A s and Af and also the universal critical amplitude ratios 
A K " I A x and A^j/A^,. The results are summarized in Table II. 

One can compare our results with the available estimates. According to a Monte Carlo 
simulation on periodic cylinders nx n x 128 with n — 4, 6, 8, and 10 (Ref. P5]), the critical 
FSS amplitude of the spin-spin inverse correlation length equals A s = 1.30(3). Our estimate 
A s = 1.26(5) is consistent with the one above. 

Information concerning the absolute amplitude of the free energy can be extracted from 
the published data by the indirect route. Indeed, in accordance with the calculations |26|| 
carried out in the quantum limit of a three-dimensional Ising model, A s /Af = 3.671(6). 
Using the above estimate A s = 1.30(3) one finds Af = 0.354(9). The value Af = 0.37(3) 
given in Table II agrees with this estimate. Note that the critical FSS free-energy amplitude 
of the three-dimensional Ising lattice with the shape of periodic cubes is equal to A C f Ube = 
0.625(5), Ref. f23| . For a comparison note also that in the two-dimensional space the free- 
energy amplitudes of periodic Ising strips and squares are equal correspondingly to n/12 = 
0.261799 and ln(2 x / 4 + 2" 1 / 2 ) = 0.639 911 (see, e.g., Ref. §). 

We do not know from the literature any estimates for the universal ratio A K » / A x in 

hh 

three dimensions. One can note that out of all the four quantities discussed in this section, 
the combination A » I A y has the best convergence in n. 

K hh A 

Finally, for the universal combination A x (a) / 'A\ which is the finite-size cumulant ratio 
^oo, there exists an estimate only in the order of magnitude: ~ 3 (Ref. f2~7[| , see also the 
reviews [[|). We have succeeded in obtaining this quantity to an accuracy of 5%. 



VI. CONCLUSIONS 

In this article, the author has presented large-scale transfer-matrix calculations for dif- 
ferent finite-size amplitudes of a fully anisotropic simple-cubic Ising lattice in the shape of 
n x n x oo bars with periodic boundaries in both transverse directions. The behavior of 
amplitude combinations (ratios) which do not contain the nonuniversal metric factors and 
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geometry prefactor was studied depending on the interaction anisotropy parameters J x /J z 
and J y /J x . It has been established that these amplitude ratios practically cease to depend 
on the anisotropy parameter J x /J z — > and, what more, this is true for wide interval, 

Jx/ Jz rS 10 

It was shown that as a function of J y j J x the amplitude ratios have a smooth extremum 
(maximum) near J y /J x = 1. As a result, the critical finite-size amplitude combinations are 
universal with respect to both anisotropy parameters in the domain < J x /J z <C 1 and 
\ Jy/Jx - 1| < 1- 

A mechanism leading to the J x / J z independence of certain amplitude ratios was pro- 
posed. By this the J x J y invariance together with the analyticity of the scaling functions 
explains the existence of a smooth extremum in the amplitude combinations at J y /J x = 1. 

In the case of fully isotropic interactions (J x = J y = J z ) for which the high accuracy 
values of critical coupling and critical free energy are available, the better estimates have 
been found for the universal critical finite-size amplitudes of the spin-spin inverse correlation 
length and singular part of the free energy per site, as well as for the universal amplitude 
ratios An /A Y and A w u\/A\. 

K hh A K ' X 
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APPENDIX: MATRIX ELEMENTS OF V 1 , V 2 , V 3 , AND V A IN THE 
QUASIDIAGONAL REPRESENTATION OF V 



The sizes of upper and lower subblocks in the expansion ( |T8| ) of TM for a cylinder 
n x n x oo are equal between themselves (Ni = N 2 ) by odd n and unequal (Ni ^ N 2 ) by 
even n. Consider the cases of even and odd n separately. 

Matrix elements of subblocks V^ 1 -* and for the bars 3 x 3 x oo and 4 x 4 x oo have 
been done in Jl5| and I. For the 2 x 2 x oo cluster, subblocks and have, respectively, 
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the sizes 5x5 and 2x2 As for the 4 x 4 x oo cluster, their matrix elements are given 
by 



and 



t/W _ a . n. . 

ij l 3 l 3 



(Ala) 



where 



T/( 2 ) _ A -rt ■ 
ij l 3 l 3 i 



exp 



11 ml + mfjK x + \{m\ + m^AT, 



(Alb) 



(A2) 



and 



G tf = gtf> + 2 £ gf j) cosh(2 S A 2 ) 

s=l 



(A3a) 



G^ = 2^#)sinh(2 S K 2 ). 



(A3b) 



The coefficients rij, m", m*, and mj for the periodic cylinder 2 x 2 x oo are 
n l = {2,8,2,2,2}, < = {4, 0, -4, -4, 4}, = {4, 0, 4, -4, -4}, m f = {4, 2, 0, 0, 0}. 



(A4) 



As in the case of 4 x 4 x oo cluster, the coefficients gyi' satisfy again the condition 



9$* + 2 £ 9P = mini 



1 2 
2 



(A5) 



s=l 



For gyfi with s ^ one has 



11) 1 21) 1 22) 1 31) 32) 1 33) 1 41) 42) 1 
43) 44) 1 51) 52) 1 53) 54) 55) 1 . 

Finally, the coefficients of a cluster 2 x 2 x oo equal 



(A6) 



11) 1 21) 1 22) 1 . (A7) 



Taking (from [DJ and I) the basis functions of irreducible representations we find the 

(12) 

matrix elements of subblocks Vf for bars nxnxoo with even n (i.e., for cylinders 2 x 2 x oo 
and 4 x 4 x oo): 

(Vi 12) h = \A l3 (mfi i3 + rmGij) , (A8) 

where rrii = by i > A^; here and below we regard, for defmiteness sake, Gij = when % or 
j > N2. Analogously one evaluates for subblocks and V^: 

(V}%j = 222!^ij [K + m ))Ga + 2m i m j Gf ij ] (A9a) 

and 

= ^M(™% + + 2m i m i G ii ] . (A9b) 

The matrix elements of subblock V3 ^ are 

(^ 3 (12) )ii = ^ ^ (3m * + m ?) G « + mi (m * + 3m ?)^'] • ( A1 °) 
Lastly, for subblock V4 we find 

(V 4 (1, )« = ^yA y [(mf + 6m?m| + + ^m^m? + m?)G y ] . (All) 

In the case of a cluster nxnxoo with odd number of layers n, the formulas are more 
uniform because the same set of coefficients enters into matrix elements for subblocks 
of both representations. Introducing the auxiliary quantities 

|(n 2 +l) 

E l^lexpps-^^sign^)] (A12a) 

and 



8=1 



4 = E |^' ) |exp[-(2 S -l)^sign(^))], (A12b) 



l(n 2 +l) 



s=l 
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we obtain for even (k) terms of transfer-matrix expansion: 

1 



k )ij 



m, t + m 3 ) k H\f + K - mj) k H^ 



(Al3a) 



and 



(Yk - 2 kk\ Aij 



[mi + rrijfHlf - (m; - m j ) k Ht ) 



i.i 



(Al3b) 



where equal ( |A2]) but, of course, with their sets of coefficients. By odd k, i.e. for subblocks 

(12) 

Vj. , the expressions for the matrix elements have the form 

1 



11 ~ 2 k k\ Al] 



mi + mjfHlJ' - (m, - mj)*H>j 



(A14) 
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TABLES 

TABLE I. Critical-point amplitude ratios A^i /A x and A x (a)A s /A x as a function of anisotropy 



parameters J x /J z and J y /J x - 


Approximat: 


l 1 j_ 

ions by cluster 


pairs (2,3) and (3,4). 




Jy/ Jx Jxj Jz 




(2, 3) pair 






(3, 4) pair 






A " /A Y 


A xW A s /A 2 x 


k B T c / J z 


A U 1 Ay 


A xW A s /A 2 x 


0+ 1.0 


2.367640 


65214 


1 92695 


2.320811 


0.48946 




0.5 


1.699790 


0.65173 


1.92388 


1.662411 


0.48991 


1.44730 


0.1 


0.921068 


0.65279 


1.92832 


0.910794 


0.49034 


1.44919 


0.01 


0.513815 


0.65309 


1.92968 


0.510586 


0.49039 


1.44940 


0.001 


0.345015 


0.65310 


1.92972 


0.343461 


0.49039 


1.44940 


0.1 1.0 


2.845802 


1.34451 


3 60933 


2.817633 


1.32110 


S 48345 


0.5 


1.979483 


1.35368 


3.62760 


1.959011 


1.34212 


3.54079 


0.1 


1.024601 


1.38006 


3.70054 


1.019932 


1.35995 


3.59108 


0.01 


0.550499 


1.38661 


3.71950 


0.549216 


1.36225 


3.59760 


0.001 


0.362466 


1.38685 


3.72021 


0.361871 


1.36233 


3.59781 


0.5 1.0 


3.819394 


1.65412 


4 55222 


3.739735 


1.67447 


4 59324 


0.5 


2.533705 


1.65621 


4.54634 


2.487136 


1.68196 


4.61381 


0.1 


1.211125 


1.67095 


4.59073 


1.199032 


1.68882 


4.63492 


0.01 


0.611620 


1.67545 


4.60591 


0.608159 


1.68976 


4.63782 


0.001 


0.390357 


1.67564 


4.60656 


0.388821 


1.68979 


4.63793 


1.0 1.0 


4.685960 


1.68710 


4.67454 


4.581044 


1.70507 


4.70713 


0.5 


3.024529 


1.68407 


4.64512 


2.960469 


1.71114 


4.72047 


0.1 


1.366300 


1.69711 


4.68102 


1.350375 


1.71743 


4.73841 


0.01 


0.658794 


1.70231 


4.69869 


0.654587 


1.71832 


4.74101 
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0.001 0.410971 



1.70256 



4.69953 



0.409173 



1.71835 



4.74110 



1.5 



1.0 

0.5 

0.1 

0.01 

0.001 



5.418190 
3.442977 
1.494583 
0.695761 
0.426622 



1.68220 
1.67334 
1.68607 
1.69307 
1.69344 



4.66797 
4.61261 
4.64298 
4.66651 
4.66777 



5.310630 
3.368233 
1.475982 
0.691216 
0.424737 



1.69188 
1.69831 
1.70726 
1.70870 
1.70875 



4.66481 
4.67598 
4.70162 
4.70596 
4.70612 



2.0 



1.0 

0.5 

0.1 

0.01 

0.001 



6.068336 
3.819394 
1.608392 
0.727256 
0.439633 



1.66945 
1.65412 
1.66577 
1.67510 
1.67563 



4.63471 
4.55222 
4.57383 
4.60471 
4.60652 



5.975643 
3.739735 
1.587593 
0.722535 
0.437724 



1.66826 
1.67447 
1.68729 
1.68970 
1.68979 



4.58676 
4.59324 
4.63016 
4.63763 
4.63792 
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TABLE II. Estimates of the universal critical amplitude combinations for the fully isotropic 
three-dimensional simple-cubic Ising lattice in the geometry of infinitely long cyclic parallelepipeds 
with the square cross section. For all quantities, with the exception of Af, the three-point extrap- 
olations were performed by the Shanks' transform. 24 In the case of Af as a realistic estimate, the 
last term of finite sequence has been taken. 



n 



A s 



Ai 



A <J A x 



A^/Al 



2 
3 
4 



1.458347 
1.353169 
1.302661 



0.451886 
0.409959 
0.378507 



1.784594 
1.763848 
1.755304 



3.497834 
3.661609 
3.759097 



oo 



1.26(5) 



0.37(3) 



1.749(6) 



3.9(2) 
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FIGURES 

FIG. 1. Amplitude ratio A k » /A x vs the anisotropy parameter J x /J z by different fixed values 
of Jy/Jx- Periodic cylinder 4 x 4 x oo at critical temperatures Tc 3 ' 4 \ 

FIG. 2. Amplitude ratios A s /Af (averaged over J x j ' J z = {1, 10 _1 , 1CT 2 , 1CT 3 }) as a function of 
Jy/Jx- Shown is are the approximations by (2,3) and (3,4) cluster pairs — dashed and solid lines, 
respectively. Open circle on the ordinate axis is a value of A s /Af(= 3) for the two-dimensional 
Ising model. 

FIG. 3. Average amplitude ratio A » j A x vs J y /J x \ periodic cylinder 4 x 4 x oo at critical 

hh 

temperatures Tc 3 ' 4 \ Full circle on the ordinate axis corresponds to the A K >/ j A x at J y /J x = + , 

hh 

while open one is a value of that ratio in the two-dimensional Ising lattice. 
FIG. 4. The same as in Fig. 3 but for the amplitude ratio A X (4)A S /A X . 
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